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For the strongly correlated topological insulator SmBe we discuss the influence of a 2 x 1 recon¬ 
struction of the (001) surface on the topological surface states. Depending on microscopic details, 
the reconstruction can be a weak or a strong perturbation to the electronic states. While the former 
leads to a weak backfolding of surface bands only, the latter can modify the surface-state dispersion 
and lead to a Lifshitz transition. We analyze the quasiparticle interference signal: while this tends 
to be weak in models for SmBe in the absence of surface reconstruction, we find that the 2x1 
reconstruction can induce novel peaks. We discuss experimental implications. 


I. INTRODUCTION 

The material SmBe has attracted significant attention 
during the last years, as it has been proposed^^^ to re¬ 
alize a three-dimensional (3D) topological Kondo insula¬ 
tor (TKI). In this fascinating class of materials, a topo¬ 
logically non-trivial bandstructure emerges at low ener¬ 
gies and temperatures as a result of Kondo screening of 
strongly correlated / electrons.^ Similar to weakly cor¬ 
related topological insulators (TIs), the surfaces of 3D 
TKIs display topological surface states with Dirac dis¬ 
persion and spin-momentum locking. 

Recent theory^~^ suggested SmBe to be also a topolog¬ 
ical crystalline insulator (TCI), its bandstructure being 
characterized by three non-zero mirror Chern numbers, 
^^^= 0 ’ ^kz='K^ defined on planes in the 3D Bril- 

louin zone (BZ) which are invariant under mirror opera¬ 
tions. This allows to further classify topological surface 
states according to their mirror-symmetry eigenvalues. 

On the experimental front, a number of results ob¬ 
tained on SmBe, e.g., from transport studies,quan¬ 
tum oscillation measurements,^^ angle-resolved photoe¬ 
mission spectroscopy (ARPES),^^^^^ and spin-resolved 
ARPES^^ appear consistent with the presence of Dirac- 
like surface states. However, doubts have been raised 
about the proper interpretation of ARPES data.^^”^^ 

An important class of experiments, typically used to 
verify the topological nature of TI surface states, em¬ 
ploy scanning tunneling spectroscopy (STS): the surface- 
state quasiparticle interference (QPI) signal allows one 
to deduce characteristic wavevectors for spin-conserving 
elastic scattering due to defects and thus directly probes 
the spin-momentum locking of surface electrons. Unfor¬ 
tunately, conclusive QPI results on SmBe are lacking to 
date. In fact, surface-sensitive probes such as STS face 
the problem that SmBe surfaces are hard to characterize 
experimentally: Non-reconstructed (001) surfaces are po¬ 
lar, such that surfaces typically reconstruct. A frequent 
case is the 2x1 reconstructed (001) surface, as observed 
in large areas in STS experiments^^’^^. Such a surface is 
non-polar, rendering it more favorable for the observation 
of topologically non-trivial surface states. This implies, 
however, that a proper theoretical modelling must take 


into account the surface reconstruction, i.e., the fact that 
one half of the top Sm rows are missing - this has not 
been done in the theory literature to date. 

It is the purpose of this paper to fill this gap: We 
shall study the effects of surface reconstruction of the 
electronic states of SmBe on the level of tight-binding 
models. In particular, we shall focus on the fate of the 
topological surface states and on signatures in STS and 
QPI experiments. Our main results can be summarized 
as follows. Not unexpectedly, the importance of the 
surface reconstruction depends strongly on the effective 
strength of the electronic reconstruction potential. The 
latter depends on various microscopic details which can¬ 
not be reliably extracted from our simplified modelling, 
such that we rely on assumptions here. In the case that 
the reconstruction potential is weak, its effects on QPI 
are minor despite the backfolding of surface bands, i.e., 
the effect of reconstruction can - for most purposes - be 
ignored. If, however, the reconstruction potential is suf¬ 
ficiently strong, a Lifshitz transition - a change in the 
Eermi-surface topology as function of the Eermi energy 
- of the surface states is expected. Its general proper¬ 
ties can be described by appealing to the presence of 
a mirror plane which prevents band hybridization along 
certain high-symmetry directions. This is in analogy to 
other TCI materials^^”^^ (where a similar scenario ap¬ 
plies to non-reconstructed surfaces). Importantly, such 
2x1 reconstruction can entirely modify the tunneling 
signal, and new QPI peaks of unidirectional character 
may appear. 

The remainder of the paper is organized as follows. We 
will start with an introductory discussion to SmBe and its 
surfaces in Section II. Section III will then qualitatively 
discuss the effect of surface reconstruction on the surface 
states. Concrete numerical results for bandstructures as 
well as tunneling and QPI spectra will be given in Section 
IV. Section V summarizes the results of our work. 


II. SmBe: GENERAL REMARKS 

To set the stage, we summarize key aspects of SmBe 
surfaces and the tight-binding modeling of both bulk and 
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surface electronic states, and we also highlight general 
aspects of surface reconstruction. 

A. SmBe (001) surface 

(001) Surface terminations. SmBe crystallizes into a 
simple cubic lattice, with lattice constant ao = 4.13A. 
Crystals cleave preferentially along the (001) direction. 
This surface, however, presents many challenges from 
both experimental and theoretical point of view, since it 
can come with many different terminations 
Be^^, Sm reconstructed^^’^^, single six or eight B 
donut-like^^, disordered^^’^^). All of these have differ¬ 
ent electronic properties, and are both difficult to con¬ 
trol experimentally and to describe in a coherent fashion 
theoretically. 

In particular, simple Sm and Be terminated surface are 
polar, which has lead to the suggestion that the observed 
surface conductance might be carried by non-topological 
surface states^^. STS experiments^^however, show 
that, while these nominally polar surface regions can be 
observed, they are somewhat rare and small, being elec¬ 
trostatically unstable. The most common situation is 
instead the one of a reconstructed surface, in which, on 
average, one half of the terminating atoms are missing, 
restoring electrostatic neutrality. 

2x1 reconstructed (001) surface. Even though this re¬ 
construction is most often disordered, large regions of or¬ 
dered, Sm terminated, 2x1 reconstructed (001) surfaces 
can be observed^^’^^. In these areas every second Sm row 
on the top layer is missing. Fig. 1. Due to the absence 
of macroscopic surface charges, this is the experimental 
scenario which most closely resembles the idealized (001) 
surface which is usually described by tight-binding-based 
theoretical models^”^’^’^^'^^. 

Impact of reconstruction on surface states. The re¬ 
moval of one half of the atoms in the top layer cannot 
be considered a harmless process; this is especially true 
if we consider surface states, which mostly live in the few 
top layers. Due to relaxation effects and the removal of 
neighboring atoms, even layers below the reconstructed 
one are expected to be influenced by the process. Here, 
the penetration depth A of surface states plays an impor¬ 
tant role: if X/oq ^ 1, these states live on many planes, 
and a strong perturbation on one or a few of them has a 
weak impact; if instead X/oq ^ 1, they live on just a few 
top layers, and a strong perturbation there can have a rel¬ 
evant effect. Using a simple model^^, it has been found 
that the penetration depth is of the order of a few atomic 
layers, suggesting the second possibility to be relevant. 

B. Modeling 

Even though it is generally agreed that the destruc¬ 
tion of one or more of the top atomic layer(s) via disor¬ 
der simply shifts topological surface states deeper in the 



Figure 1. Pictorial view of the (a) non-reconstructed (001) 
surface and (b) of the 2x1 reconstructed one; we only show 
Sm atoms, and in (b) we highlight with two different colors 
the two sublattices a and b. 


bulk^^’^^, not much is known about a strong but periodic 
perturbation, as the 2x1 reconstruction that we consider 
here. It is therefore the aim of this paper to discuss the 
effect of this scenario on the electronic surface states. To 
this end, we shall employ an eight-orbital tight-binding 
model, following Refs. 3, 6, and 29 (see the Appendix for 
details). 

Bulk mean-field approach. We recall that SmBe is a 
strongly interacting material, in which local / moments 
are coherently screened by conduction d electrons to lead 
to a Kondo insulating phase. At low temperatures, this 
physics can be captured using a single-particle descrip¬ 
tion with renormalized parameters'^: in particular, the 
/ kinetic energy is strongly suppressed, by a factor that 
dynamical mean-field theory (DMFT) and the Gutzwiller 
approximation find to be ^ Simpler theories 

such as mean-field slave bosons^^ can qualitatively re¬ 
produce this reduction, but with a smaller renormaliza¬ 
tion factor ^ 2^’^^. Instead of describing interaction- 
induced renormalizations explicitly, here we will work di¬ 
rectly with renormalized parameters'^ that are designed 
to match qualitatively the computed bandstructure and 
the measured ARPES dispersion. 

Role of surfaces and Kondo breakdown. In presence of 
spatial inhomogeneities, the interaction-induced renor¬ 
malizations will be position-dependent^^. In particu¬ 
lar, Kondo screening can be appreciably reduced on a 
surface, leading to an effective decrease of the Kondo 
temperature. This can generate a so-called “Kondo- 
breakdown” , in which surface moments are no more 
Kondo screened: they can be free, order magnetically, or 
be possibly still screened but with a Kondo temperature 
smaller than the bulk one. Decoupled moments lead to 
an increase of the Dirac velocity, and to a displacement of 
the Dirac energy deep in the valence band, which resem¬ 
bles experimental results indicating “light” rather than 
“heavy” surface states. In Ref. 36 the Kondo tempera¬ 
ture was estimated to drop by roughly a factor 3, as a 
consequence of surface Sm atoms having one first-nearest 
neighbor (NN) Sm atom less than in the bulk. In pres¬ 
ence of a 2 X 1 reconstruction atoms on the top row will 
have three neighbors less: in this case, the reduction of 
the hybridization, hence of the Kondo temperature, at 
the surface would be even larger. 
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Atom 

1st NN 

2nd NN 

3rd NN 

A(bulk) 

6 

12 

8 

H(non-rec. surface) 

5 

8 

4 

C(rec., top row) 

3 

4 

4 

D(rec., below the missing row) 

5 

10 

8 

F(rec., below the existing row) 

6 

10 

4 


Table L Number of Sm neighbors for Sm atoms (A) in the 
bulk, (B) on the (001) non-reconstructed surface, (C) in the 
hrst [(D,E): second] row on the 2x1 reconstructed (001) 
surface; different positions A-E are illustrated in Fig. 2(a-d). 


However, an analysis^’^^ of tight-binding models built 
to fit ab-initio calculations indicates that the bulk gap is 
mostly generated by 2nd NN hybridization, and in Ref. 
37 a model with 3rd NN hybridization only seemed to 
reproduce well neutron-scattering data. This then leads 
to a more complex scenario, where the reduction of the 
hybridization of the surface atoms is less drastic. We 
summarize the count of different NNs in Table I. While 
we will not model Kondo breakdown explicitly, we will 
return to this aspect in the next paragraph. 

Additional surface scattering potential. To formally 
eliminate half of the atoms on layer 1, we introduce a vir¬ 
tually infinite reconstruction potential Vj-ec acting every 
second row (in practice we take Vrec = 100 eV). With Vrec 
being active, we need to introduce an additional surface 
scattering potential AE = AE^ = AEf of a few hundred 
meV, acting on the d and / orbitals of the remaining 
atoms in layer 1, in order to qualitatively reproduce the 
surface-state dispersion as seen by ARPES (large Dirac 
velocity, Dirac energy in the valence band). For / elec¬ 
trons, which have a bandwidth of a few meV, this value 
turns out to be very close to the infinite scattering po¬ 
tential of the Kondo-breakdown scenario AEf = oc, in 
which surface / electrons are no more part of the coher¬ 
ent Fermi liquid. Even though quantitative details might 
vary between our approach and the Kondo breakdown 
one, most of the general features we are going to describe 
are general, and independent of microscopic details. We 
also remark that, while a large AEf is needed to repro¬ 
duce light surface states, AE^ is introduced to further 
adjust the position of the Dirac energy, hence the Fermi 
momentum kp (that we try to obtain as close as possible 
to experiments), and may not be needed for some choices 
of the parameters (see Section IV C). 


C. Qualitative effects of reconstruction 

In the following we choose a coordinate system as in 
Fig. 2, with z perpendicular to the surface and y being 
parallel to the surface rows. 

Surface Brillouin zone. A 2 x I reconstruction doubles 
the size of the primitive surface unit cell in x direction 



Figure 2. Sketch of the Sm sites for a (001) surface (a)-(c) 
without reconstruction and (b)-(d) with a 2 x 1 reconstruc¬ 
tion. In the first case, of C 4 v symmetry, four mirror planes are 
present (corresponding to operators M^, My, Mx-y, Mx+y), 
while in the second case, of C 2 v symmetry, only Mx and My 
remain; in this second case we depict the two sublattices 
with different colors. Corresponding two-dimensional Bril¬ 
louin zones (e) without reconstruction (“large” BZ) and (f)- 
(g) with reconstruction; in this second case the BZ is shrunk 
along the x direction (“small” BZ), and the X cone is folded 
at r (f); however, for weak reconstruction, it is useful to use 
the original large BZ (g), with weak band replicas appear¬ 
ing as a consequence of backfolding (dotted ellipses). In the 
non-reconstructed case high-symmetry points are labeled as 
r = (0,0), V = (7r,0), X' = (0,7r) and M = (7r,7r), and 
mirror planes are kx = 0, = 0, kx — n, ky — tt, ky — H-fc, 

ky = —kx, which are pairwise equivalent by C^v symmetry. 
In the reconstructed case, high-symmetry points are labeled 
as r = (0,0), X 2 = (7r/2,0), Y = (0,7r) and M 2 = (7r/2,7r), 
and mirror planes are kx — t), ky — 0, kx — ky = tt . 


and introduces two sublattices, a and b with and without 
top row, respectively. Fig. 2(a-d). The two-dimensional 
(2D) surface BZ, originally a square for the (001) sur¬ 
face, shrinks by one half in x direction according to the 
reconstruction wavevector Q = (^,0), Fig. 2(e-f). Here 
we will use the term “2nd BZ” for the difference between 
the original (“large”) and the folded (“small”) BZ. 

Spectral weights. A (weak) reconstruction causes a 
backfolding of bands which leads to replicas (or shadow 
bands) shifted by Q w.r.t. the original bands and with 
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reduced spectral weight. Experimentally, a weak replica 
of the X cone has been observed at ARPES, 

even though the origin of this signal is not totally clear 
(for example, having polycrystalline samples, one would 
expect to observe the replica of the X' cone at E, too). 

Within our model calculations, we can quantitatively 
define spectral weights in the following way. The wave- 
function '0^(k, z, a, m) of a single-particle state with en¬ 
ergy En for a system in slab geometry depends on the mo¬ 
mentum k = {kx, ky) in the small BZ (—7r/2 < kx < 7r/2, 
—IT < ky < ir), on the layer index z, on the orbital in¬ 
dex Of, and on the sublattice index m = 0(a), 1(6). Its 
spectral weight in the large BZ is given by: 


w(n^ K, z, a) = 


Kxm 


m=0,l 


'0n(k,2;,a,m) 


( 1 ) 


with K = {Kx^Ky) from the large BZ (—tt < Kx^Ky < 
tt), and k its backfolded partner, k = K mod Q. We 
can define a global (i.e. orbital-, layer-, and momentum- 
integrated) energy-resolved weight separately for the 1st 
or 2nd BZ: 

w^^\E)= w{n,K^^\z,a)L{E-E„), 

,z,a 

w^^\E) = Y. ( 3 ) 

n,K 2 rid^ 2 ;,Ck: 

where L is a Lorentzian kernel, and means \Kx\ < 
7r/2, while means \Kx\ > 7r/2. We can also define 

a layer-resolved weight 


'^z{E)= y] w;(n,K,2;,Q!)L(£; - (4) 

n,K,Ck: 


and decompose it into contribution from the 1st and 2nd 
BZ by limiting the sum over K accordingly, to yield 
and w^^\E), with wf^\E) + w^^^{E) = 
w,{E), and Ez^z{E) = p{E) = w^^\E) E w^^\E), 
where p{E) = L{E — E^) is the density of states. 

Effective reconstruction potential To quantify the ef¬ 
fect of reconstruction on surface states, we introduce the 
qualitative notion of “effective reconstruction potential” 
(ERP) strength which is based on spectral weights. ERP 
is a dimensionless quantity between 0, implying no effect, 
and 1, implying maximal effect. 

We start from an idealized situation, in which the re¬ 
construction potential only affects wavefunctions in the 
first layer z = 1. On this layer, the weight in the 1st and 
2nd BZ are equal, w{n^ z = l^a) = w{n, z = 

l,(a) = w{n,'K,z = l,(a)/2, as only the m = 0 term in 
Eq. (1) survives. In all other layers, the weight (of a 
single state) will be totally either in 1st or in the 2nd 
Brillouin zone, according to where it was without recon¬ 
struction. Hence, the ratio of the layer-integrated weights 
in the 1st and 2nd BZ (i.e. the relative weight of replicas) 
can be used to infer the weight of surface states on the top 
layer. Specifically, for a band located in the 2nd BZ with¬ 
out reconstruction (like the X cone), with w{n^ z > 


I, a) = 0, > l,a) = u;(n,K,2: > l,a) 0, 

the normalized weight on the first layer can be written 
as: 


w 


/ _ 
1 — 


Wi 

P 


yyBZl UiBZ2 ’ 


0 < rei < 1. 


( 5 ) 


We note that the above assumption is not strictly valid: 
We found that the weight of the X cone can be non- 
negligible in the 1st BZ even for deeper layers z > 1, even 
though no explicit reconstruction potential is present 
there, see Section IV. Nevertheless, we shall use w[ from 
Eq. (5), evaluated for a surface state at the Eermi level, 
as a measure of the ERR We will see in Section IV that 
this value can also be cone-dependent. 

Row strong is the reconstruction? We now briefly dis¬ 
cuss whether existing experimental or theoretical data 
are sufficient to infer the effective strength of the recon¬ 
struction potential. 

Experimentally, ARPES data indicate that band repli¬ 
cas are weak and hence the ERP, as measured by Eq. (5), 
is small^^’^^. This is, however, difficult to reconcile with 
our model calculations which, under essentially all cir¬ 
cumstances, show that the removal of every second Sm 
atoms in the surface layer is a strong perturbation to the 
surface states, with large ERP (see Section IV for de¬ 
tails). Possible reasons for this discrepancy are: (i) only 
a fraction of the surface probed by ARPES features an or¬ 
dered 2x1 surface reconstruction, thus reducing the ERP 
in the collected ARPES signal, or (ii) surface states in the 
real material tend to be expelled from the first layer(s) 
by some other mechanism, such as Kondo breakdown, or 
(iii) the penetration length of surface states is too small in 
the model calculations. Issue (i) may be investigated with 
small-spot ARPES experiments, and we will not discuss 
this further. Eor (ii), we found that the introduction of 
the surface scattering potential, XE^ and AEj, to model 
this expulsion often enhances rather than decreases the 
ERP, but this is a model-dependent statement. We note 
that ARPES may not probe the full weight of surface 
states if their penetration length is larger than that of the 
ARPES probe, but this is unlikely to explain the above 
discrepancy. The issue (iii) of correctly modelling the 
surface-state penetration is delicate: The problem of the 
excessive reduction of the bulk gap using renormalized 
parameters (with respect to plain DET ones, noted for 
example in Ref. 38) has a direct impact here: if one wants 
to match the experimental ARPES gap ^ 20 meV using 
renormalized parameters, one has to increase the value of 
the hybridization, effectively decreasing the penetration 
length of surface states, and, as a consequence, enhancing 
the strength of the ERP. The problem is further compli¬ 
cated by Ty states^’^^”^^, whose inclusion would reduce 
the value of the bulk gap (in our model, we will ignore 
these states, but a more detailed analysis would need to 
include them). 

We conclude that we are not in the position to fully 
resolve this discrepancy, and further work - both theo¬ 
retical and experimental - is needed to do so. Instead, we 
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will take a pragmatic point of view and analyze the fate 
of the surface states provided that the ERP is sufficiently 
strong to generate observable effects. 


III. SURFACE STATES WITH 
RECONSTRUCTION 

In this Section we qualitatively analyze the effect of 
surface reconstruction on the in-gap surface states, based 
mainly on symmetry arguments. Numerical examples, 
taking into account details of microscopic modelling and 
testing these predictions, will be provided in Section IV. 



A. Symmetries and dispersion 

Energies elose to the Dirae point. If the Dirac energies 
were close to the Fermi level, the Fermi wavevector kp 
of the cones were small, and essentially the only effect 
of reconstruction would be the folding of bands into the 
small BZ, as shown in Fig. 2(f)-(g). (The same applies 
at all energies for weak ERP.) 

Overlap of X eone with its repliea. Away from the 
Dirac energies, Dirac cones and their replicas will start 
to overlap with each other. ARPES data on SmBe show 
that the only relevant overlap is the one of the X cone 
with its replica: 2kF of the X cone along the X — V 
direction is estimated to be^^ around 0.80 A“^, while the 
original BZ has a halfwidth 7r/ao = 0.76 A“^. Thus, 
we expect that, due to the reconstruction, the X cone 
should overlap and hybridize with its replica for energies 
close to the Fermi energy; we will ignore similar effects 
which might affect the X' (now Y) and F cones at higher 
energies. 

Hybridization, role of mirror symmetry, and Lifshitz 
transition. We now discuss general aspects of the cou¬ 
pling between the X cone and its replica which are largely 
independent of microscopic details; see Fig. 3. A key in¬ 
gredient is the reflection symmetry with respect to the 
xz plane, described by operator My: States with ky = 0 
are invariant under My (Ref. 5). Along this direction 
the two cones have the smallest distance, and the over¬ 
lapping states have opposite mirror eigenvalues ± i (their 
exact value depends on the sign of the mirror Chern num¬ 
ber but is irrelevant for the argument). Hence, the 

cones cannot hybridize for ky = 0 (as long as mirror sym¬ 
metry is preserved), instead they cross at E = Fcross- 

In contrast, hybridization is allowed for ky ^ 0. We 
observe that for energies E\^ corresponding to momenta 
k < 7 r/ 2 , the X cone deforms {E = Ei in Fig. 3) 
along the ky = 0 direction, until a Lifshitz transition 
at E = Emerge < Ecross Is achicvcd: the X cone and 
its replica merge at ( 7 r/ 2 ,/c^), where in general 7 ^ 0 . 
As a consequence, a pocket centered at X 2 = (' 7 r/ 2 , 0 ) 
is created; we note that the minimum along X 2 - M 2 is 
actually a saddle point. When the two cones cross at 
E = Ecross along the ky = 0 direction, the X2 pocket 


Figure 3. Lifshitz transition for the X cone, (a) Dispersion 
along the f — A 2 — M 2 lines (the f cone is a spectator). 
At E = Emerge the X conc merges at (kx,ky) = (7r/2,±/c°) 
with its own replica, centered at F, giving rise to the Lifshitz 
transition, while at E = Ecross the two cones cross along the 
ky = 0 line: along this direction they cannot hybridize, having 
different mirror-symmetry My eigenvalues. We also report 
mirror eigenvalues ± i, and draw the SEV with green arrows 
(for solid arrows the relation to mirror eigenvalues is univocal, 
for dotted ones is model dependent), (b) Schematic evolution 
of the X cone and of its replica across the transition; with a 
dotted ellipse we show how the cones would evolve without 
reconstruction potential Free- 

area vanishes. At higher energies E > Ecross, the X 2 
pocket grows again, such that a new Dirac cone is cre¬ 
ated at the X 2 pocket from the crossing of the two orig¬ 
inal Dirac cones. We note that for vanishing scattering 
potential, the two cones always cross without hybridiz¬ 
ing: this yields Emerge = Ecross, and the two branches 
along the X 2 — M 2 direction (kx = 7r/2) are degener¬ 
ate. Finally, for the case with reconstruction and broken 
mirror symmetry, a gap is opened at X 2 . 

Analogy to TCIs. The hybridization behavior in 
Fig. 3 is similar to what happens on the (001) sur¬ 
face of face-centered-cubic TCIs such as Pbi_a,Sna,Se^^, 
Pbi_a,Sna,Te^^, and SnTe^^. For these materials, for 
which the standard Z 2 indices are zero and which are 
not TIs in the usual sense, mirror Chern numbers are 
non-zero, and predict two pairs of Dirac cones along the 
kx = 0 and ky = 0 directions of the 2D BZ. These cones 
are at (tt ± (5/^, 0) and (0, tt ± (5/^), where the BZ is defined 
by —IT < kx^ky < IT. Each pair of cones merge around 
(tt, 0) or (0, tt) through the same mechanism that we have 
described here, and which has been also predicted^ for 
the (110) surface of SmBe. The similarity lies in the 
non-trivial topology of the band structure and the pres¬ 
ence of mirror planes, which create states with different 
mirror eigenvalues which cannot hybridize along specific 
directions. 

However, in the TCI case, no surface reconstruction is 
needed to generate the pair of cones. Also, the role of 
symmetries is reversed: in standard TCIs mirror symme¬ 
try protects cones at (0,7r ± Sk), while parity invariants 
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predict two Dirac cones at (0, tt). In our case we identify 
(0,7r) with X 2 and (0,7r ± Sk) with F and X, but now 
mirror symmetry protects the cones at X 2 , while parity 
invariants predict cones at F and X. 

ARPES experiments on TCIs have verified^^’^"^’^^ the 
surface-state behavior as depicted in Fig. 3. While a sim¬ 
ilar verification might be possible in SmBe, we note a few 
caveats: Visible hybridization effects require the ERP to 
be strong; as noted in the previous Section its magnitude 
is unknown. Moreover, the features might be beyond the 
present ARPES energy and momentum resolution, and 
finally the cone crossing might occur above the Fermi 
energy, rendering ARPES ineffective. 


B. Direction of spin 

So far, the analysis of crossing Dirac cones was based 
on the presence of a mirror plane, but was otherwise in¬ 
dependent of microscopic details. We now take a step 
further and compute mirror-symmetry eigenvalues and 
the expectation value of the electronic spin operator. 
As shown in Refs. 6, 7, 40, and 41, mirror-symmetry 
eigenvalues are determined by the sign of mirror Chern 
numbers (see Appendix). While four TCI phases with 
^ 1=0 = ±2, = +1, = ±1 are possible in 

principle, we restrict the analysis to the ~ 

= 7 r ~ =k ~ phase, which - as we have 

argued in Ref. 6 - is the one of experimental relevance 
for SmBe. In Fig. 4(a) we report the mirror-symmetry 
eigenvalues for the non-reconstracted surface.^ 

Mirror-symmetry eigenvalues with reeonstruetion. In 
the reconstructed case, there are four distinct mirror 
planes in the 2D BZ, namely = 0,7r/2 and ky = 0,7r; 
the former commute with M^,, the latter with My (see 
Fig. 2(b)-(d)-(f)). We can therefore assign mirror eigen¬ 
values to states crossing these planes. For planes = 0, 
ky = 0^ ky = TT everything is equivalent to the non¬ 
reconst r acted case, and we can transfer the eigenvalues 
directly from Fig. 4(a) into Fig. 4(b). For energies below 
Emerge: ^o statcs cross kx = 7r/2, and no further work is 
needed. Above Emerge: however, four states are crossing 
the kx = 7r/2 plane, two for ky > 0, and two for ky < 0. 
Concentrating on ky > 0 we find that below Ecross both 
states have eigenvalue — i, see Fig. 4(c), simply because 
they belong to the same band, see Fig. 3(a). In contrast, 
above Ecross: Fig. 4(d), ouc of the ky > 0 states has eigen¬ 
value + i, and the other — i, see Fig. 3(a), since they now 
belong to different bands; the same happens for ky < 0, 
see Fig. 4(c,d). 

Spin expeetation value (SEV). Since mirror-symmetry 
eigenvalues are not directly measurable, we now concen¬ 
trate on the SEV, which depend on these eigenvalues 
and is accessible in SP-ARPES experiments. For atomic 
d and Fg states, that we consider here, the relation be¬ 
tween mirror eigenvalues and SEV is straightforward, and 
it is the same one as for free spins^: if the eigenvalue is 





Figure 4. Mirror-symmetry eigenvalues ± i and SEV for (a) 
the non-reconstructed and (b-d) the reconstructed BZ at dif¬ 
ferent energies, see Fig. 3. Mirror eigenvalues and SEV not 
reported in panels (c,d) are the same as in panel (b). 


— i, the SEV points along -\-x for Mx and -\-y for My^ and 
viceversa for eigenvalues + i. 

This holds also for lattice states with a definite mo¬ 
mentum, when symmetry operations act like in the non- 
reconstructed case. This applies to states kx = 0, ky = 0, 
OT ky = TT - hence, on these planes, the SEV is identi¬ 
cal to the non-reconstructed case - but not for states at 
kx = 7r/2 because this is not a BZ mirror plane in the 
absence of reconstruction. 

Let us analyze the latter case of kx = 7r/2 in more 
detail. For the two sublattices a and b we define mirror 
eigenstates according to Mx\a~^) = +i|a+), Mx\b'^) = 
-l-i|6+), and we build Bloch states |k+) and |k^) with 
k = (7r/2, /cy). When we act with the mirror operator Mx 
on |k+) we get Ma,|k+) = -hi |k+), but for states on sub¬ 
lattice b we pick up a minus sign: Ma,|k^) = — i |k^) [the 
opposite would happen if we took a mirror operator Mx 
centered on b rather than on a in Fig.2(b-d)]. This can 
be realized by observing that Mx^ when centered on sub¬ 
lattice a, sends states a in the n primitive cell into states 
a in cell (—n), while b states at n are sent into b states 
at (—n — 1): so, the effect of Mx on jk^) states carries 
an additional factor, which is —1 for kx = 7r/2. 

As a consequence, states |k) with Ma,|k) = +i|k) and 
k = (7r/2, /c^) are composed by states |a+) on the sublat¬ 
tice a, and by states \b~) on sublattice 6, so there is no 
general rule for their SEV given the mirror eigenvalues. 
We can however state that, if we restrict the analysis 
to the top layer, where no b states are present, eigen¬ 
values — i corresponds to a SEV pointing along , and 
viceversa for states with -hi. When including all layers, 
however, we have to resort to numerical diagonalization. 
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In our model we find that, close to the Lifshitz transition, 
the total contribution from other layers is opposite to the 
one from layer 1 , so the SEV for — i states points along 
—x; at higher energies, however, we observe a transition 
to a situation where the SEV for — i states points along 
(see Eig. 3(a)), in agreement with the results from 
layer 1 . 

IV. NUMERICAL EXAMPLES AND QPI 

In this section we provide numerical examples, using 
explicit tight-binding calculations, of the scenario de¬ 
picted in the previous Section, together with simulated 
QPI patterns. 

A. Introduction to QPI 

Generalities. QPI can be used as a probe of the topo¬ 
logical character of surface states^^”^^, owing to their spin 
texture: non-magnetic impurities cannot induce transi¬ 
tions between states which are Kramers partners. Also, 
when mirror symmetries are present, impurities which 
do not break these symmetries cannot induce transi¬ 
tions between states with opposite mirror eigenvalues^^. 
In the following, we will exclusively consider impurities 
which do not break time-reversal and mirror symmetries; 
more complicated impurities may induce additional QPI 
peaks^^’^^ which we will not discuss here. 

In general, peaks in the QPI signal are expected for 
pairs of stationary states^^ (states which have the same 
tangent in the BZ to the isoenergy contour), for which 
no selection rule forbids an elastic scattering process. Eor 
TI surface states, no peaks are expected from intracone 
scattering as long as the dispersion is linear, since pairs 
of stationary points are time-reversal conjugate in the 
Dirac Hamiltonian^^. Peaks can arise when the disper¬ 
sion deviates from linear, as it is well known to happen 
in TIs such as Bi 2 Ti 3 due to hexagonal warping^^. 

Intereone seattering and QPI on SmB^ non- 
reeonstrueted (001) surfaees. Intercone scattering 
is less restricted by selection rules and may hence 
give rise to strong QPI peaks^^. However, for the 
non-reconstructed (001) SmBe surface, we showed in 
Refs. 6 and 29 that the experimental spin pattern, corre- 
spending to the = + 2 , = + 1 , C+= -1 

mirror Chern numbers and to a positive winding num¬ 
ber w = sgn(C^ ^^) on the X cone, leads to a 
strong suppression of the intercone scattering, implying 
virtually no observable QPI peaks. This most likely 
explains why STS experiments on SmBe have not found 
noticeable QPI signals to date^^’^^. 

QPI in SmB^ with reeonstruetion. A key question is 
whether QPI peaks can appear as a consequence of sur¬ 
face reconstruction. As noted above, the question makes 
sense only for sufficiently strong ERP, otherwise the non- 
reconstructed picture, which shows no peaks, still holds 


to a good approximation. We can expect new QPI scat¬ 
tering channels in the energy region around the Lifshitz 
transition: first, because for energies below the transi¬ 
tion an appreciable warping of the X cone is achieved, 
second, because above the transition the topology of the 
Eermi surface changes, with the appearance of the new 
pocket at X 2 . Concrete calculations are required to see 
whether these prerequisites yield a noticeable QPI signal 
- those will be presented in what follows. 

B. Numerical results: Significant warping 

We employ the model described in the Appendix, taken 
from Refs. 6 and 29, which is designed to match as closely 
as possible the calculated bulk bandstructure, and to re¬ 
produce the experimentally observed surface-state dis¬ 
persion and spin pattern. Without surface reconstruc¬ 
tion it yields no sharp QPI peaks, and consequently the 
appearance any QPI peaks can be traced to reconstruc¬ 
tion effects. Eor a specific set of parameters which pro¬ 
duces significant warping near the cone crossing, we re¬ 
port the band dispersion for a slab geometry in Eig. 5, 
the corresponding density of states (DOS) in Eig. 6 , and 
constant-energy cuts through the ARPES and QPI sig¬ 
nals in Eig. 7. Here, we define the ARPES signal as the 
sum over the weights from all orbitals^^, and add the 
signal from layers 1 and 2 , 

A{E, K) = X E (6) 

n,Q: 2;=1,2 

expressed using the weight defined in Eq. (1). 

Dispersion and DOS. In the dispersion, Eig. 5, we can 
see two cones at T, and one cone at Y. Eor completeness, 
we present in Eig. 5(a) the case without reconstruction, 
Kec = O 5 and in Eig. 5(b) the case with reconstruction 
but no surface scattering potential, AE = 0. In what 
follows we mainly concentrate on Eig. 5(c), which shows 
the case with both reconstruction and surface scattering 
potential. 

The weight distribution in Eigs. 7(a-d) shows that the 
Y cone is weakly affected by the reconstruction, and lives 
mostly in the 1st BZ: for this cone the ERP is weak. 
In contrast, the two cones present at T are strongly hy¬ 
bridized with each other, and have appreciable weight 
both in the 1st and in the 2nd BZ even when consid¬ 
ering the signal from deeper layers: for these cones the 
ERP is strong. Their origin is from T and X in the non- 
reconstructed BZ; in analogy to the experimental situa¬ 
tion we identify as “T” the cone with Dirac energy close 
to the top of the valence band and smaller Eermi momen¬ 
tum, and as “V” the cone with lower Dirac energy and 
larger Eermi momentum (see Eig. 3). 

We can observe that the X cone undergoes the Lifshitz 
transition that we described in the previous section, while 
the T cone remains inert. The saddle point at Emerge 
gives rise to a peak in the DOS (strictly speaking, a log¬ 
arithmic divergence), as shown in Eig. 6 . Parenthetically, 



15 20 



Figure 5. Dispersion of a 15-layer slab along the rX 2 M 2 Yr 
path for (a) Kec = 0, AE = 0; (b) Kec = 100 eV, AE = 0; 
(c) Vrec = 100 eV, AE = —0.5eV. The Fermi energy has been 
placed in the middle of the bulk gap. Line colors encode the 
added weight on layers 1 and 2; for case (c) we also label 
surface states with mirror-symmetry eigenvalues ± i. We can 
observe two cones at F and one at T. In (c), Emerge ~ 4meV 
and F^cross ~ 7 meV, the latter corresponding to the Dirac- 
cone crossing at W- Horizontal dashed lines correspond to 
the energies of the cuts in Fig. 7. 


we note that no cone crossing is observed inside the gap 
for significantly smaller values of the surface scattering 
potential AE (see for example Fig. 5(b)). 

QPL To illustrate QPI effects, we consider a single 
impurity (corresponding to the dilute limit) and employ 
the standard T matrix technique to compute the Fourier- 
transformed local DOS in the presence of elastic scatter¬ 
ing; for details see the Appendix. We report two sets of 
results: In Fig. 7(e-h) we show the QPI signal from layer 

1, which has equal weight in the 1st and in the 2nd BZ 
due to the missing b atoms. In contrast, in Fig. 7(i-l) 
the signal is constructed from the sum of the local DOS 
in layers 1 and 2, the latter scaled down by a factor 10, 
to account for its larger distance from the tip. A more 
refined analysis might include interference contributions 
from d and / channels^^’"^^ as well as from layers 1 and 

2, but would lead to the same qualitative results. 

It can be observed that far from the Lifshitz transi- 


—I-1-1- 

1 St layer d (atom C) 
1st layer f (atom C) 
2nd layer d (atom D) 
2nd layer f (atom D) 



Figure 6. Local DOS for d and / orbitals on atoms C and 
D from Fig. 2, for parameters as in Fig. 5(c). The peak at 
Emerge ~ 4 mcV (arrow) originates from the Lifshitz transi¬ 
tion; the peak around — 17meV is the bulk / band. Dotted 
vertical lines around ±9 meV denote the bulk band gap. The 
/ DOS on atom D is rescaled down by a factor 200. 


tion no QPI signals, apart from an incoherent peak at F, 
are present (Fig. 7(a)); this is the same result as without 
reconstruction. At higher energies, slightly below the 
Lifshitz transition a weak intracone scattering signal is 
generated as a consequence of cone warping (wavevector 
qi in Fig. 7(b)), while above the transition two scattering 
channels appear to be active, both involving the newly 
formed X 2 cone. In the first one the scattering is towards 
the newly formed band coming from the merging of the X 
cone with its replica, denoted q 2 in Fig. 7(d), and which 
involves states with the same mirror eigenvalue with re- 



Figure 7. (a-d) ARPES signal from layers 1 and 2, Eq. (6), 
(e-h) QPI signal from layer 1, and (i-1) QPI signal from layers 
1 and 2 (divided by 10), all for parameters as in Eig. 5(c). Ar¬ 
rows illustrate scattering wavevectors. Different columns cor¬ 
respond to different energies: (a,e,i): far below the Lifshitz 
transition, yielding no coherent QPI signal; (b,f,j): slightly 
below the transition, where a QPI signal at qi arises from 
warping (the solid and the dotted green arrows are equiva¬ 
lent up to an inverse lattice vector); (c,g,k): at the transition; 
(d,h,l): above the transition, with new scattering channels 
q 2,3 involving the new pocket at X 2 . The QPI signal is nor¬ 
malized to the LDOS at the Eermi energy in the hrst layer. 
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E=-3meV 
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E=+4meV 



E=+7meV 
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Figure 8. Same as Fig. 7(a-h), but for set of parameters with 
a reduced ratio between hybridization and / kinetic energy, 
see text. This case yield no structured QPI signal. 


spect to Mx on the = 7r/2 plane; this is in analogy 
with the results of Ref. 45 for Pbi-a^Sn^^Te. The second 
one is towards the Y cone (qs); in some cases we also 
observe a weak signal coming from the scattering to the 
F cone. We stress, however, that the intensities (both 
absolute and relative) of these new QPI peaks depend on 
the exact choice of model parameters, as well as on the 
parameters used to simulate the tunneling process from 
the tip to the surface. We also note that these peaks 
arise primarily from the spectral weight in the first layer, 
which is mostly of d character. Since d states contribute 
a small fraction of the total weight of the surface states, 
their QPI peaks may - depending on details of the tun¬ 
neling process - be masked by an incoherent background 
arising from / states. This is illustrated in Fig. 7(i-l) 
where we have included a contribution from the second 
layer (which does carry / weight); as a result, most QPI 
peaks dramatically lose contrast. Even under these cir¬ 
cumstances it is still possible to observe a weak signal 
from the qi channel. 


C. Numerical results: Other cases 

We have repeated the same calculations for various 
other sets of parameters, and we describe a few repre¬ 
sentative cases in what follows. 

Kondo breakdown. To model Kondo breakdown in the 
top layer, we take AEf = 100 eV, much larger than all 
other energy scales. Combined with the same AE^ = 
—0.5eV, we obtained numerical data (not shown) that 
are almost indistinguishable from the example shown be¬ 
fore in Fig. 7. We remark that a finite AE^ is needed 
in our set of parameters to increase kp and allow for the 
Lifshitz transition, but is not be needed for all sets of 
parameters (see below). 

Reduced hybridization vs. kinetic energy. If we reduce 
the hybridization with respect to the / kinetic energy, 
while keeping approximatively the same gap (t^ ^ 

V 0.5i;), the QPI signal strength is strongly reduced. 
Fig. 8. In particular, for energies below the Lifshitz tran¬ 
sition, we observe no sign of nesting. This can be reason- 


Figure 9. Same as Fig. 7(a-h), but for a weak reconstruction 
potential Wc = 0.2 eV, AE = —0.4eV. Weak replicas of the 
cones can be observed in the isoenergy contour (a-d), and no 
appreciable QPI signal is present (e-h), even though a signal 
corresponding to q 2 in Fig. 7 starts to be visible. 



Figure 10. Same as Fig. 7(a-h), but for for a moderate recon¬ 
struction potential Wc = 0.5 eV, AE = —O.SeV. 


ably linked to the slightly increased penetration length 
of surface states, which makes reconstruction a weaker 
perturbation. In general, the appearance of QPI peaks 
and in particular their strength is found to be parameter- 
dependent. 

Surface scattering potential AE^ = 0. For a similar 
set of parameters (see Appendix) we are able to get a 
kp close to the experimental value and realize a Lifshitz 
transition with AEf = 100 eV and AEd = 0. The result¬ 
ing QPI signal (not shown) is similar to Fig. 8. 

Small reconstruction potential. For smaller reconstruc¬ 
tion potential F^ec = 0.2 eV and FFec = 0.5 eV, Figs. 
9 and 10, the ARPES signal displays weak replicas of 
the cones appearing translated at the reconstruction 
wavevector Q = (tt, 0). In this case the QPI signal is very 
weak and almost indistinguishable from the case F^ec = 0? 
which, as remarked, shows no peaks. 


V. CONCLUSIONS 

In this paper we have studied the role of a periodic 
2x1 surface reconstruction on the topological surface 
states of SmBe using a simple tight-binding model. We 
qualitatively distinguish two cases, according to whether 
the effective reconstruction potential (ERP) acting on 
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the surface states is weak or strong. Weak ERP only 
produces a backfolding of surface bands, such that weak 
replicas of the original Dirac cones appear shifted by the 
reconstruction wavevector Q = (tt, 0). In contrast, strong 
ERP induces a particular crossing of Dirac cones, accom¬ 
panied by a Lifshitz transition of in-gap states and the 
formation of a new Dirac cone, protected by mirror sym¬ 
metry, at the edge of the small BZ. We have provided a 
numerical example for the case of strong ERP, and shown 
that new QPI peaks can appear as a consequence of this 
transition. 

The ERP itself depends mostly on the weight of sur¬ 
face states on the first atomic layer(s), which in turn 
depends on microscopic details such as the penetration 
depth A of surface states, or the presence of a surface 
scattering potential. Our tight-binding approach cannot 
fully predict the strength of the ERP because detailed 
information on the structure of the reconstructed surface 
would be required. Ab-initio calculations could help, but 
to our knowledge no systematic studies have been con¬ 
ducted on reconstructed SmBe surfaces. Experimentally, 
existing ARPES data - showing weak replicas only - sug¬ 
gest a weak ERP which, however, is difficult to reconcile 
with our analysis: Using reasonable model assumptions, 
we tend to generically obtain strong ERP. Possibly, the 
ERP is stronger than what appears from ARPES exper¬ 
iments, for example if the probed surface area contains 
large non-reconstructed (or disordered) regions which in 
turn do not contribute to band backfolding. We recall 
that some other correlated materials with strong peri¬ 
odic modulations also fail to display strong bandfolding 
effects in ARPES, presumably due to quenched disorder, 
one example being Lai. 8 -xEuo. 2 Sra,Cu 04 ^^. 

Based on our results, we suggest that small-spot 
ARPES on 2 X 1 reconstructed surfaces as well as care¬ 
ful QPI studies, searching for additional peaks appearing 
above the Lifshitz transition, could clarify and further 
elucidate the surface-reconstruction effects in SmBe- 
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Appendix A: Tight-binding modeling 
1. Bulk model 


To model SmBe, we only retain Sm atoms, and we 
build a tight-binding model out of the da, 2_^2 - 
quadruplet in the d shell, and of the Pg quadruplet 

in the / shell, where = y^l ± |) + y^l ^ D’ 

jPg^^T) = I ± kinetic energy terms 

for d states {tf = ((i ^2 |iL|(i^ 2 )ooi, = (d ^2 |iL|(i^ 2 )iio, 

where the suffix 001 or 110 denotes the direction along 
which the matrix element is considered), two kinetic 
energy terms for / states {t{ = (Pg^^ + |iL|rg^^+)ooi, 
tl = (Pg^^ + |iL|rg^^-h)iio), two hybridization terms 

(^;i = (4^ t |ff|rf+)ooi, V 2 = t |ff|rf4)iio), 
and one on-site energy difference between d and / states 
Cd — ef. The kinetic energy is fixed to give at the three 
X = (0, 0, tt) = (0, TT, 0) = (tt, 0, 0) points a minimum 
in the d shell, with da, 2_^2 symmetry at (0,0,tt), and 
a maximum in the / shell; this leads to band inver¬ 
sion at the X points, and to the topological Z 2 indices 
(z/q, 5 ^ 3 ) = (I 51515 !)• The hybridization is chosen 

to be in the = +2, = +1, = -1 

phase, where are mirror Chern numbers^”^: 


2 ^ r 

^ E E 

a,6=1 n=l^^ 


(f]<.{daU+{'k)\dbU+{'k)), 


(Al) 

where M|iz+(k)) = +i|iz+(k)) and k lying in the plane 
BZ which is invariant under the symmetry operator M 
{M=Mz when BZ is = 0 or = tt, M = M^-y 
when BZ is kx = ky). This leads to a positive winding 
number on X cones^ as observed experimentally^^, and 
is relevant for the QPI signal and the spin structure of 
surface states, but not for their dispersion. 

The Hamiltonian 


H — Hkin T Hdyhr 


is composed by the kinetic energy term: 


Hkin — 

and by the hybridization: 


Hd 0 
0 Hf 


Hhybr — 


0 

-iU/ 


0 )' 


(A3) 


(A4) 


The kinetic energy Hi, i = d/f, is diagonal in the 
(pseudo)spin index, with basis dx 2 -y 2 and ^^2 when 

i = d, and Pg^^ Pg^^ when i = f; the hybridization 
Hdf is nondiagonal in the (pseudo)spin index, with basis 

dx^-y^ dx 2 -y 2 dz 2 t ^nd dz 2 ^ for rows, and Pg^^T, 
Tg^^-, for columns. 

Their explicit form is as follows {cx = coskx^ Cy = 
cosky, Cz = cos kz, Sx = sin/ca,, Sy = sin ky, Sz = sin/c^, 
5-[- Sx i i Sy, CS^ ^y^x T i CxSy^. 
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^ + ^{c^ + Cy)(t\/2 + Czt\) - Cy) [-t\/2 +c^t\) 

y V^{Cx — Cy) {—t\l2 + 02 ^ 2 ) + (c® + Cy) (ti/2 + 02 ^ 2 ) + '2Czt\ + 4CxCyt2 


i = d/f, 


(A5) 


^ (^x “t“ ^y)^z 

Zs+{vi/2 + CzV2) 
V^(Cx Cy)^2^2 

V3s_(-vi/2 + c^V2) 


3s_(vi/2 + €2^2) 
-3v2(Cx + Cy)Sz 
V3s+{-Vi/2 + CzV 2 ) 

-V3(Cx - Cy)SzV2 


Parameters are taken from the tight-binding calcula¬ 
tions of Refs. 28 and 29, with a slave-boson-like renor¬ 
malization of the / kinetic energy by and of the hy¬ 
bridization by 6, with ^ 0.08, and fine tuned to better 
match ARPES experimental data. For Figures 5, 6 and 
7 we increased the value of the hybridization by roughly 
50% with respect to this estimation to better highlight 
QPI features; this however does not affect appreciably 
the ERP, which would be in any case large. In partic¬ 
ular we employ: = 1.8 eV, tf = —0.76eV, = 

0.2eV, t{ = 3.2 meV, ^2 = —l.OmeV, vi = —84meV, 
V 2 = 24meV, which lead to a bulk gap A = 18meV, to 
= —lOmeV with respect to the chemical potential, 
and to a level occupation rid = 0.51, rif = 3.49, which 
corresponds to a Sm^-^+ valence; the chemical potential 
is set in the middle of the bulk gap. In Fig. 11 we show 
the dispersion for a slab without reconstruction; the same 
data is shown in Fig. 5(a) in the small BZ. 


> 

CD 

E 

LU 



0.5 

0.4 

0.3 

0.2 

0.1 

0 


Figure 11. Bandstructure of the SmBe model employed in 
the paper, in a 15-layer slab geometry without surface recon¬ 
struction, plotted along the f AMf path in the large surface 
BZ. 


V3s-{-vi/2 + CzV2) 
Sz[2vi + {Cx + Cy)v2] 
s+{vi/2 + CzV 2 ) + 4v2CS+ 


\/3s+(-t>i/2 + CzV 2 ) 

-\/3{Cx - Cy)SzV2 

s_(vi/2 + CzV 2 ) + 4i;2CS_ 

— Sz[2vi + {Cx + Cy)V2] 


\ 

(.A6) 

/ 


2. Surface reconstruction 

We model reconstruction by applying a large scatter¬ 
ing potential Idee = 100 eV on sites every second row of 
the top layer; we neglect any relaxation effect that the re¬ 
moval of these top atoms can cause. We also apply a sur¬ 
face scattering potential AE = AEd = AEf = —0.5eV, 
which is set to align the energies of the Dirac cones to 
approximatively reproduce a situation close to the exper¬ 
imental one, allowing for the Lifshitz transition of surface 
states. The penetration length A at the Fermi energy is 
about two atomic layers for all cones. 

In the Kondo-breakdown scenario we set AEf = 
100 eV for / states, but still retain a finite AEd = 
—0.5eV for d states; in this case no appreciable modi¬ 
fication of the dispersion or of the QPI signal is observed 
with respect to Figures 5 and 7 of the main text. 

In Fig. 8 we use t{ = 6meV, = —3meV, vi = 
—42 meV, V 2 = 12 meV, AE = 0.1 eV, and keep the other 
parameters fixed. 

We are able to achieve the Lifshitz transition with 
AEd = 0? when t{ = 4.8 meV, = —1.8meV, vi = 
—46 meV, V 2 = 13 meV, and the other parameters fixed. 
For this parameter set we do not show results, but QPI 
patterns are very similar to those of Fig. 8. 


3. QPI 

QPI figures are generated with a 400 x 400 mesh on a 
15-layer slab, using the scattering matrix technique (see 
Ref. 35 for technical details) with an artificial broaden¬ 
ing of (5 = ImeV (Figs. 9 and 10: 5 = 2meV). We take 
a 30meV scatterer in both d and the / channels in the 
first two layers in sublattice a, but results were found to 
be qualitatively independent of the choice of the scatter¬ 
ing potential. As remarked earlier, quantitative results 
depend on the exact choice of the parameters'"^. 

In the QPI figures, we only report the Fourier trans¬ 
form of the impurity-induced piece of the LDOS, i.e., 
structural peaks at momenta (0, 0) and (tt, 0 ) (the recon¬ 
struction wavevector) are not shown. 
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